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ABSTRACT 


Electrical  network  analysis  from  associated  system 
matrices  generally  results  in  determinants  which  are 
awkward  to  handle  because  computation  becomes  labori¬ 
ous  with  an  appreciable  number  of  meshes.  This  paper 
is  a  study  of  several  ladder  -type  networks  where  recur  - 
sion  formulas  for  the  system  determinants  are  solved 
by  the  method  of  finite  differences.  It  appears  that  a 
broad  class  of  networks  are  amenable  to  this  type  of 
analysis  and  certain  generalizations  of  the  subject 
method  are  set  forth. 


PROBLEM  STATUS 


This  is  an  interim  report  on  this  problem;  work  is 
continuing. 


iv 


A  FINITE  DIFFERENCE  SOLUTION  OF  RECURRENT  NETWORKS 


INTRODUCTION 

It  is  well  known  that  electrical  networks  can  be  analyzed  from  the  standpoint  of  their 
associated  system  matrices.  In  general,  however,  the  resulting  determinants  become 
awkward  to  handle.  Should  a  circuit  have  an  appreciable  number  of  meshes,  the  computa¬ 
tions  become  extremely  laborious.  The  present  paper  is  a  study  of  several  ladder- type 
networks,  uniform  and  otherwise  where  recursion  formulas  for  the  system  determinants 
are  obtained  which  are  then  solved  by  the  method  of  finite  differences.  It  appears  that  a 
broad  class  of  networks  are  amenable  to  this  analysis.  With  the  ladder-type  circuits  as 
points  of  departure,  certain  generalizations  of  the  method  are  set  forth. 


THEORY 

A  basic  recurrent  ladder  network  is  illustrated  in  Figure  1. 

Z,  Z,  Z,  Z, 


By  definition, 


Figure  1 


Zx  +  ,?Z2  =  S 
and  Z2  =  B 


Then  the  loop  equations  become,  in  matrix  notation, 


(1) 

(2) 


(3) 


1 


2 
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where  designating  the  first  matrix  on  the  right  hand  side  as  Jm.  J , 

M..  =  S 
11 

M.  .  =  M.  .  .  =  -B 

i,  1  +  1  i+  1,  l 

Mj.  =  0  (j  +  i,  j*i+l) 

Equation  (3)  could  also  be  written 

M  ■  M  0 

the  argument  denoting  the  number  of  meshes  in  the  system. 

With  S  +  0,  let 

1  =0 

The  case  of  S  =  0  is  considered  later. 

It  follows,  therefore,  that 

[e]  =  S  [Q(n)]  [i] 


where 


Q(n) 


k-A 

-A 

0 


-P 


-P 

1 


(4) 

(5) 

(6) 

(7) 


(8) 


(9) 


(10) 


A  complete  analysis  of  the  circuit,  Figure  1,  requires  evaluation  of  the  determinant 
of  [M(n)j  ,  its  zeros,  and  its  first  cofactors. 


Defining  the  determinant  of 
and  the  determinant  of 
then 


M(n)  =  D(n) 

Q(n)  =  A(n) 

D(n)  =  Sn  A(n) 


(ID 

(12) 

(13) 


SOLUTION  OF  A  (n) 

Consider  the  determinant  A  (n). 

Expansion  by  the  first  row  (column)  results  in 

A(n)  =  A(n-l)  -  f?  A(n-2)  (14) 

Where  A  (n-1),  A  (n-2)  are  determinants  of  the  same  form  as  A  (n)  but  of  dimensions 
(n-1),  (n-2),  respectively. 


NAVAL  RESEARCH  LABORATORY 


3 


Assume  a  solution 

A(n)  =  CAn  (15) 

where  C  and  A  are  constants. 

Direct  substitution  into  (14)  yields 

An  An-1  2An-2 

A  =  A  -p  A  (16) 


or 


yielding 


A2  -  A  +  p2  =0 

A  =  11  V1‘4PX 
2 


(17) 


(18) 


Defining  Ai  and  A2  as 

A  -  UlSmi  a,  =  lliSZ 

2  2 

The  complete  solution  is  then 

A(n)  =  Ct  Axn  +  CaA2n 

since  (14)  is  a  linear  finite  difference  equation  and  (16)  through  (18)  justify  choice  of 
bases  AL  and  A2.  Where  Cx  and  C  2  are  arbitrary  constants,  it  should  be  noted  that 

Aj  +  A2  =  1  (211 

A,  A2  =  p*  (22) 


(19) 


(20) 


To  evaluate  the  constants  Ct  and  C2,  two  boundary  conditions  are  required. 
By  definition  of  A  (n) 


A(l)  =  1 
A(2)  =  l-pa 


Using  (20)  through  (24),  it  follows  that 


-  A  A2 

C.  =  T-i 4-  and  C  = 


<i  A2-At  - 2  A2-At 

Hence  the  complete  solution  becomes 


a<“>  ■  "VA  f 


A.n+1-A 


,n+1] 


(23) 

(24) 


(25) 


(26) 
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Solving  (37),  and  using  the  boundary  conditions  (38)  and  (39), 

D(n)  =  -jn  z£_  [(-l)n+1-l] 
2 


Hence 


D(n)  =  0  if  S  =  0  and  n  is  odd,  and 

D(n)  =  (-l)n/^  Z  n  when  S  =  0  and  n  is  even. 
2 


(40) 

(41) 

(42) 


ZEROS  OF  D(n) 

To  obtain  the  natural  modes  of  vibration  of  the  system,  which  are  required  in  analyzing 
the  transient  behavior,  the  zeros  of  D(n)  must  be  evaluated. 

If  B  =  Z2  =0  then  from  (3)  D(n)  =  Sn  =  Z* 

Thus,  unless  ZL  and  Z2  are  zero  simultaneously,  Z2  =  0  does  not  yield  a  zero  of  D(n). 

The  preceding  analysis  shows  that  only  zeros  of  S*n*?  ^+1)0 nee(j  be  considered. 

smhp 

If  sinh  (n+l)0  =  0,<f>  =  j^-  . 

However,  since  sin  kn  =  0,  k  =  0  and  k  =  n+1  must  be  excluded. 

Hence  zeros  occur  when  k  =  1,  2,  ...,  n 

Specifically,  then,  since  cosh  j  9  =  cos  9,  zeros  are  determined  by  p  =  l/2  sec-^— 

^(43) 

It  should  be  noted  that  equation  (41)  yields  a  possible  source  of  zeros,  namely  if  n  is 
odd  and  S  =  0. 

This  case  corresponds  to  choosing  k  =  n-t-l/2  in  (43),  and  solving  for  p.  The  solution  is 
p  =  oo,  i.e.  S  =  0. 

The  circuit  of  Figure  2  may  be  considered  as  an  example. 


R  R  R 


Figure  2 

Z i  =  Ra,Z2  =  l/£p  where  p  is  the  required  angular  frequency  mode. 


s 
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„  ,  k tt  ,  RCp  2  r  kff 

Zeros  occur  when  cos  =  1  +  — jp- ,  or  p  =  cos  ^  -1 

GENERAL  STEADY  STATE  BEHAVIOR 


To  determine  the  steady  state  currents,  it  is  necessary  to  evaluate  the  first  cofactors 
of  D(n).  If  the  driving  voltage  is  placed  in  the  first  mesh,  the  relevant  cofactors  are  Di  .(n) 
(j  =  1,  2,  ...,  n)  and  for  applied  voltages  in  the  k-th  mesh,  Dkj(n)  ( j  =1,  2,  ...,  n).  1 


By  a  Laplacian  expansion  about  the  first  (r-1)  rows  of  Dr  r+k 

Dr,r+k  (n)  =  Bk°(r-1)  D(n-r-k) 

with  the  convention  that  D(0)  =  1 
Since  D(n)  is  symmetric 


(n),  it  is  apparent  that 
(44) 


Dkj(n)  =  Djk(n) 


(45) 


0r  Dr+k  r  ^  =  ^  D^r_1)  (D-n-r-k) 


(46) 


Hence,  for  a  voltage  et  applied  to  the  first  mesh,  the  current  in  the  k-th  mesh  is 

given  by  D  (n) 

,  kw  ^k-l  D(n-k)  . 

\i  D(n)  ei  B  D(n)  i  (47) 


and  substituting  the  appropriate  values  from  (2)  and  (35), 


ii  = 

sinh  (n-k+1)  <!>  0 

Z2  sinh  (n+l)<£  1 

(48) 

Similarly  for  voltage  in  the  jth  mesh 

Bj_k  D(k-l)  D(n-j) 

(k<j) 

(49) 

ikf 

D(n) 

e. 

1 

or 

ikj  = 

sinh  k  (f>  sinh  (n-j)<£ 

Z2  sinh  (n+1)  0  sinh^ 

e. 

1 

(k<j) 

(50) 

and 

iki  = 

Bk-j  D(j-l)  D(n-k) 

e. 

(k>j) 

(51) 

KJ 

D(n) 

3 

which  is  equivalent  to 

V 

sinh  j  <t>  sinh  (n-k)<^ 
Z2  sinh  (n+l)<A  sinh# 

®j 

(k>j) 

(52) 

Finally,  for  voltages  ev  e2,  ...en  in  all  the  loops 


*k 


_  1 

Z2  sinh  (n+1)  <f>  sinh  <f> 


'k 

2 

j-1 


e.  sinhj^sinh  (n-k)0  + 
1  t 


n 

£  e,  sinh  k  <p  sinh  (n-t)  4> 
=  k+1  1 

(53) 
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GENERAL  TRANSIENT  BEHAVIOR 

The  preceding  analysis  carries  over  to  the  transient  case  as  well.  For  the  steady 
state,  Z2  and  <t>  are  functions  of  the  impressed  angular  velocities,  a>.  for  the  transient  case 
they  are  the  same  functions  of  the  natural  modes  p..  The  amplitudes  of  the  transient  cur¬ 
rents  are  proportional  to  the  cofactors  Dj^  (n).  If  j  ^  k,  the  cases  of  j>k  and  j<k  must  be 
distinguished  as  in  equations  (50)  and  (51). 


TAPERED  LADDER  STRUCTURE 

The  preceding  analysis  applies  to  a  uniform  ladder  structure.  In  this  section  the 
applicability  of  the  method  to  a  nonuniform  structure  will  be  indicated. 

Consider  the  network  of  Figure  3. 


Z,  KZ,  K2Z,  K^'z, 


Let  Sj  be  the  self  impedance  of  the  j-th  loop 


Sx=Zl+  (1+k)  Z2 


S.  =  Kj_1  S, 
i  . 


The  system  determinant,  Dy  (n),  becomes 


(54) 

(55) 


(56) 


(57) 
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and  expanding  by  the  first  row  of  the  determinant  results  in 

D'  (n)  =  kn^2  [Sl  D'  (n-1)  -  k  z\  D'  (n-2)j  (58) 

which  can  be  solved  by  the  same  method  as  that  for  the  uniform  ladder  structure. 


Let 

Then 


Where 


p'  =  Z^A 

D'  (n)  =  kn/2  Z,n  Q' 


(n) 


(59) 

(60) 

(61) 


Expanding  Q'  (n)  by  its  first  row  (or  column) 

Q'  (n)  =  Q'  (n-1)-  kp'2  Q'  (n-2) 

Hence,  p'  j/lTreplaces  p  of  the  uniform  ladder  structure. 

The  solution  of  (62)  becomes 

DMn)»k”//'2  zf  -SinfijMlit' 
2  smh^' 


(62) 


(63) 


where 


cosh  <£' 


2  P 


X_  1  fz, 

'VZ-f£  |_2Zz 


(l+k)Zz 


2Z, 


(64) 


For  k  =  1,  the  solution  reduces  to  that  of  the  uniform  ladder.  The  currents  are  obtained 
in  a  similar  manner. 


CONCLUSIONS 

The  preceding  analysis  applies  to  networks  actuated  by  ideal  generators.  The  case 
for  generators  and  loads  of  arbitrary  impedance  can  be  readily  obtained  for  the  steady 
state;  the  transient  case  involves  the  solution  of  a  transcendental  equation  which  cannot, 
in  general,  be  expressed  in  a  closed  form  but  requires  design  curves.  With  the  more 
common  types  of  terminations,  however,  transient  solutions  are  obtainable. 
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